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Abstract 

In this paper we present algorithms for computing the topology of planar and space 
rational curves defined by a parametrization. The algorithms given here work di- 
rectly with the parametrization of the curve, and do not require to compute or use 
the implicit equation of the curve (in the case of planar curves) or of any projection 
(in the case of space curves). Moreover, these algorithms have been implemented in 
Maple; the examples considered and the timings obtained show good performance 
skills. 



1 Introduction 



The topology of planar algebraic curves, implicitly given, is a well-studied 
problem (see ([H), (fiol ). (11), (12), (13), and the more recent works (Q), (17), 



among others); more recently, the problem for space algebraic curves has also 
received certain attention (see (0), (0), Isl)). In all these works it is assumed 
that the curve is given by means of implicit equations, and the considered 
algorithms deal with the curves in this form. However, in this paper we address 
the problem, apparently not discussed up to now, of computing the topology of 
a rational curve (i.e. constructing a planar or space graph describing the shape 
of the curve) starting directly from its parametrization, without computing 
or making use of the implicit equation of the curve. This question may be 
of special interest in the field of computer-aided geometric design (CAGD), 
where many of the curves used are rational and even directly provided in 
parametric form (e.g. Bezier curves, B-splines, NURBS). 
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Perhaps the reason for the absence of previous studies in this direction is the 
common behef that if the parametric equations of a curve are available, the 
curve is easy to visualize. This is essentially true, but if the goal is to get 
a global idea of how the curve is like, then there are still some difficulties. 
On the one hand, one should previously compute a parameter interval such 
that the plotting of the curve over the interval shows the main features of 
the curve; this includes handling the case of possible missing points/branches 
(which happens if some point of the curve is generated when the parameter 
of the curve tends to infinity, see for example (0)). On the other hand, the 
plotting, as pointed out by Gonzalez- Vega and Necula in the introduction to 



(I12I). may not always provide a clear idea of the topology of the curve, and 



hence auxiliary tools for describing the shape of the curve may be of help. 

In this paper we address both planar and space rational curves. As in other 
topology algorithms, we require the input curves to satisfy certain conditions 
that can be achieved with generality. In the planar case it is required that 
the curve has neither vertical asymptotes nor vertical components, and that 
the parametrization is proper (see Section [2]). Initially, the algorithm works 
in a similar way to existing algorithms, i.e. first one computes the critical 
points and the points of the curve lying on the lines x = ai containing some 
critical point, and then one appropriately connects these points. However, the 
connection phase is carried out not in the usual way, but taking advantage of 
the fact that a parametrization is available (see Theorem [8] in Section [3]). More 
precisely, the algorithm computes how the parameter values are connected; 
so, two points are joined whenever the algorithm detects that the parameter 
values giving rise to them need to be connected. In particular, and unlike many 
classical algorithms, this strategy does not require the curves to be in generic 
position (as defined in ffl^). 

The method is specially profitable in the case of space rational curves. Existing 
implicit algorithms compute the topology of the curve by projecting it onto a 
plane (the x|/-plane, in our case), and then lifting to space the topology of this 
projection. In ([2I), this lifting phase is carried out in general by using a second, 
auxiliary projection; however, in (0), Isl) no auxiliary projection is needed. In 
any case, the lifting of the singularities of the projection is a delicate opera- 
tion. In our case, we use a similar strategy for 3D curves. However, here the 
lifting operation (which is performed without auxiliary projections) presents 
no difficulties since the space points are identified by the parameter values 
giving rise to them (previously computed when addressing the projection). In 
the case of 3D curves, our requirements are: (i) the curve has no asymptotes 
or components normal to the xy-plane; (ii) the projection onto the xy-plane 
fulfills the requirements of the 2D algorithm. 

We have implemented the algorithms in Maple 13; outputs and timings of 
several examples are given in Section O In our implementation we give to the 
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user the option of computing isolated points of the curve or not. The reason 
for this is that isolated points correspond to points generated by complex, 
non-real, values of the parameter, and therefore they may not be of interest 
for certain users; moreover^he number of isolated points is certified by means 
of Hermite's method (see (p)) and therefore it may be time-consuming. 

The structure of the paper is the following. In Section [2] we provide the neces- 
sary background on rational curves; hence, notions like properness, normality, 
critical and singular points are reviewed here, jointly with related results. In 
Section [3] we provide the algorithm for the 2D case. In Section IH the algo- 
rithm for the 3D case is given. Finally, in Section [5] we describe some details 
of the implementation, and we provide the outputs and timings of different 
examples in 2D and 3D. The parametrizations used in the examples are given 
in Appendix I and Appendix II. 



2 Background on Rational Curves 



In this section we briefly recall the background on affine rational curves that 
we need in order to develop our results. So, in the sequel we will consider an 
affine rational curve C defined by a rational parametrization 



^{t) = (Xi(t),X2(t), . . .,Xn(t)) 



' Pl{t) P2{t) Pnjty 



where gcd(pi,gi) = gcd(p2,g2) = • • ■ = gcd{pn,qn) = 1 and Pi{t),qi{t) G Z[t] 
for all i = 1, . . . ,n. In our case n = 2 or n = 3; so, we will usually write 
X, y, z instead of X2, x^. Moreover, since the parametrization is assumed to 
be real, we have that C is a real curve (i.e. that it consists of infinitely many 
real points), although for theoretical reasons when necessary we will see the 
curve embedded in C". Nevertheless, our goal will always be the description 
of the shape of its real part. 

A point Pq G is reached by the parametrization if{t) if there exists ^ C 
such that v'(^o) = Po] in this case, we will also say that to generates Pq. Notice 
that the value of the parameter generating a real point may be either real or 
complex, and that there may be points generated by several (real or complex) 
values of the parameter. In this sense, we will say that the parametrization 
ip{t) is proper if almost all points of C are reached by just one value of the 
parameter t, i.e. if (p{t) is injective for almost all the points of C. So, if (f{t) 
is proper then there are just finitely many points of C generated by several 
different values of the parameter, corresponding to the self-intersections of 
the curve. In order to check whether ip{t) is proper, we will use the following 
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criterion. Let 

Gi{t,s) =pi(t)gi(s) -pi{s)qi{t) 
G2{t,s) =p2{t)q2{s) -P2{s)q2{t) 

Gn{t, S) = Pn(t)qn{s) - Pn{s)qn{t) 
(5(^,5) =gcd(Gi,G'2,...,G„) 

Then, the following theor em, directly deducible from Proposition 7 in (jisl ) 



see also Theorem 4.30 in (jl9l ). for the planar case), holds. 



Theorem 1 The parametrization (f(t) is proper ijf G{t, s) = t — s. 

On the other hand, we will say that (f{t) is normal if every point in C is reached 
by at least one value of the parameter, i.e. if v?(C) = C. If ip{t) is not normal, 
then (see Proposition 4.2 in (Q)) there is just one point of C non- reached by 
the parametrization, namely the point 

Poo = limt-,±ooip{t) 

Notice that Poo exists if and only if deg(pj) < deg(gj) for all i e {1, 2, . . . , n}. 
Furthermore, if Poo exists, it may still be reached by some (real or complex) 
value of the parameter. If we denote Poo = (ai, a2, . . . , On), Poo is reached iff 

deg (gcd(aigi(t) - pi{t), a2q2{t) - P2{t), anqnit) - p„,(t))) > 1 

Also, observe that if Poo exists then it is obtained as the limit of a sequence of 
real points of C, and therefore it cannot be isolated. Hence, if Poo is reached 
by some value ta then it is a self-intersection of the curve, because it is a 
crossing of two branches of the curve, one corresponding to t ^ ±oo and 
the other corresponding to ta- On the other hand, if Poo exists but it is not 
reached, one can reparametrize the curve so that it is reached (see Theorem 
7.30 in fflil)). However, reparametrizations may complicate the equations of 
the curve, or bring other difficulties, like improperness or the introduction of 
algebraic numbers. Hence, in our case whenever we meet this phenomenon, 
we will understand that this reparametrization has not been performed. 

If every point of C is reachable via ip{t) by real values of the parameter one 
says that ip{t) is M-normal. We refer to (jl), (jiol ) for a thorough study of 
this phenomenon. If if{t) is not M-normal, then there exist real points P & C 
reachable only by complex values of the parameter. Moreover, the following 
result (see Proposition 4.2 in (j3)) clarifies the nature of these points. 

Proposition 2 Let (f{t) be a proper parametrization of C. Then P ^ P^o, 
P E C is non-reached by any real value of the parameter if and only if it 
is a real isolated point of C. 
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2. 1 Critical Points of Planar Rational Curves 



In the rest of the section we assume that n = 2, i.e. that C C is a real 
rational curve parametrized by ip{t) = {x{t),y(t)). Let / G be its 

implicit equation; then we have the following classical definitions: 



Definition 3 A point P E C is called: (a) a critical point if f{P) = ^{P) = 0; 

(b) a singular point, if it is critical and (-P) =0; (c) a ramification point if it 
is critical, but non-singular; (d) a regular point if it is not critical. 



One may easily see that ramification points correspond to those points sat- 
isfying that x'(t) = but y'{t) ^ 0, and that singular points correspond to 
either points where x'it) = y'{t) = 0, or to self-intersections of the curve. 
Singularities of a rational parametrization can be computed directly from the 
parametrization, without converting to implicit form. More precisely, the fol- 
lowing result holds (see Theorem 10 and Theorem 11 in (fl^). Here, we denote 
= Gi/G, G2 = G2/G, and we write M{t) = Res,(Gi,G2). 



Theorem 4 Let ip[t) be a parametrization of C, and let Pq E C be an affine 
singularity of C, reacheable by some value to E C of the parameter. Then, 
M{to) = 0. 



Remark 1 If Poo is reached by some to E C (in that case it is a self-intersection 
of the curve, and therefore a singularity, as we observed before), then to rnust 
be a root of M{t) (see Theorem 10 in lis))- 



Whenever ip{t) is proper, one may deduce that M{t) is not identically 0; 
therefore, in that situation M{t) has finitely many roots and from Theorem 
m the t-values generating reachable singularities are among these roots. Now 
let us denote the numerator of x'{t) by N{t), and let us write the square- 
free part of M{t) ■ N(t) as rh(t); also, let q(t) = lcm(gi,g2), and let m(t) = 
m{t)/ gcd(m(t), Then, the following corollary on the real critical points 
of C can be deduced. 



Corollary 5 The real critical points of C are included in the (finite) set con- 
sisting of: (i) Poo (if it exists); (ii) the real points generated by (real or complex) 
roots ofm{t). 
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3 Computation of the Graph Associated with a Planar Curve 



Let C C be a planar algebraic curve, parametrized by 

Pijt) P2{t) 
,?iW'92(t) 



ip{t) = {x{t),y{t)) = (pM^P^) , gcd(pi(t),gi(t)) =gcd(p2(t),g2(t)) = 1 



In this section we address the problem of algorithmically computing a graph 
Q homeomorphic to the curve C. In order to do so, we will follow the usual 



strategy widely used in the implicit case (see (|8|), (jl2l ). (jl3l ). (llTl)): 



(1) Compute the critical points of C (see Definition [3] in Subsection 12. ip . Let 
ai < ■ ■ ■ < be the x-coordinates of the critical points of C; also, let 
ao = -oo, am+i = +00. 

(2) Compute the points of C lying on the vertical lines x = ai, i = l,...,m 
passing through the critical points; we will refer to these lines as critical 
lines. 

(3) For i = l,...,m — 1, compute the points of C lying on the vertical line 
X = {ai + aj+i)/2; similarly for x = ai — 1, x = + 1. We will refer to 
these lines as "non-critical" lines. 

(4) Connect, by means of segments, the points of C lying on each non-critical 
line, with the points in the critical lines immediately on its right and on its 
left, respectively. 

In our case, we will take advantage of the fact that a parametrization of the 
curve is available; this will be specially useful in order to carry out step (4). 
Moreover, in order to apply the method presented in this section, we need 
that certain hypotheses are fulfilled by C. These hypotheses are introduced in 
Subsection 13. 1[ Then, in Subsection 13.21 and Subsection 13.31 we show how to 
compute the vertices and edges, respectively, of the planar graph. Finally, in 
the last subsection we provide the full algorithm. The reader may find several 
examples of the output of this algorithm in Section [5l 



3. 1 Hypotheses 



In the rest of the section, we assume that the following hypotheses are fulfilled: 

(i) (fit) is proper. 

(ii) C has no vertical asymptotes; in particular, it is not a vertical line. 

The first hypothesis guarantees that C is traced just once when following the 
parametrization <^(t). In order to check this hypothesis. Theorem [1] can be ap- 
plied. Moreover, if this hypothesis does not hold, one can always reparametrize 
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the curve (see Chapter 6.1 in (jl9l)) so that it is fulfilled. In order to check the 
second hypothesis, one can use the following result, which is easy to prove. 



Lemma 6 C has a vertical asymptote iff one of the following conditions oc- 
curs: (a) q2(t) has some real root which is not a real root ofqi{t); (h) deg(p2) > 
degfe) hut deg(pi) < deg(gi). 

If C has some vertical asymptote, one proceeds in the following way: 

• If C has no horizontal asymptotes (which can be checked by appropriately 
adapting Lemma [6]), then by interchanging the axes x and y the condition 
is fulfilled. Notice that this is an affine transformation, which therefore does 
not change the topology of the curve. 

• If C has also horizontal asymptotes, then almost all changes of coordinates 
of the type {x = X + uY,y = Y}, with /i G Q, set the curve properly 
(see Proposition 3.2 in (llll)). Observe that if (p{t) is proper, the curve 



obtained by applying such a transformation is properly parametrized by 
^^it) = ix{t)-ixyit),yit)). 



3.2 Vertices of the Graph. 



The notable points of C are the real critical points. Now from Corollary [51 we 
have that these are among the following points: 

(i) Poo (if it exists). 

(ii) The points of C generated (via ^{t)) by the real roots of the polynomial 
m{t) in Corollary O 

(iii) The real points of C generated (via ^{t)) by complex roots of m{t). 

The computation of Poo is described in Section |2l Moreover, once the real 
roots of m{t) are computed, the points in (ii) are obtained by evaluating 
x{t),y(t) at these roots. Now we consider as vertices of the graph Q not only 
these points, but also the points of C lying on the vertical lines containing the 
points in (i) and (ii). In order to compute these points, we recall the definition 
of the polynomials Gi,G2,G, introduced in Section [2l and we consider the 
polynomials 

G.(M)^£iM, G.(M)-SiM. 

^ ^ G{t,s) ' ^ ' ^ G(t,s) 

Then, given a point P^ = {xr,yr) = '^(fr), tr G M, the real roots of Gi{t,tr) 
provide the t- values of the points lying in the line x = x^; then, the coordi- 
nates of those points can be obtained by evaluating x{t),y{t) at these t- values. 
Observe that we get not only the coordinates, but also the t- values generating 
the points, via (p{t). This is important for the connection phase. 
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So, let us consider the points in (iii). If a point in (iii) is also generated by a 
real value of the parameter, then it will have already been computed as a point 
in (ii). So, if this is not the case, by Proposition [2] it is an isolated point. Now 
these points might be computed by seeking complex roots of m(t) giving rise 
(when evaluating x{t), y{t)) to real points of C. However, in the sequel we will 
provide an alternative way for carrying out this computation, that allows to 
certify the existence or non-existence of this kind of points. For this purpose, 
we denote a complex value of the parameter t = u + iv, where i"^ = —1 and 
M,f G M, and we represent the complex modulus as | ■ |. Also, we write 



Pi{u + iv) ■ qi{u + iv) 



and 



P2(u + iv) ■ q2{u + iv) 



\(l2? 1^2^ 



(a(u, v) + ih{u, v)) 



{c{u, v) + id{u, v)) 



Then the following result, that can be easily verified, holds. Here, we denote 
the result of substituting t = u + iv in q{t) = lcm(gi, ^2)5 as v). 

Lemma 7 Let Pq & C (IM."^. Then, Pq is generated by a complex value of the 
parameter to = uq + ivQ if and only if there exists Wo & M. satisfying that 
{uo,vo,wo) is a real solution of the system 

b{u,v) = 

d{u,v) = (1) 
V ■ \q{u, v)\'^ ■ w — 1 = 



In order to certify the number of real solutions of System (JT]) we apply Her- 
mite's method (see for example ({g])). However, these solutions include the 
complex values of the parameter generating real points that are also reached 
by real values of the parameter. In order to identify the existence of those so- 
lutions, we compute, also by Hermite's method, the number of real solutions 
of the system obtained by adding the following equations to System ([T]): 



\qi{u,v)\^ 
, , hiu, v) 

yit) = I L 
m{u,v)\- 

V ■ q(t) ■ w — 1 = 



(2) 



So, real isolated points of C correspond to solutions of System ([T]) which are 
not solutions of System ([2]) . 
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3. 3 Edges of the Graph. 



In this section, we address the problem of connecting the vertices of Q (to 
compute the edges of the graph). For this purpose, the idea is to introduce 
between two consecutive critical lines an intermediate "non-critical" line, and 
to connect the points of C on each "non-critical" line with the points of C on 
the critical line immediately on its right or on its left. In order to do this, we 
take advantage of the fact that a parametrization of the curve is available, 
and we connect the points just by comparing the parameters generating the 
points in the two vertical lines (one of them critical, and the other one "non- 
critical"). The idea is made precise in the following theorem. Here, we will 
consider P^o as "generated" by both +00 and —00, besides other real values 
that may also generate it; as usual, —00 (resp. +00) is considered less (resp. 
greater) than any other real number compared with it, and —00 < +00. This 
result is illustrated by Figure 1. 

Theorem 8 Let Xa,Xh G M satisfying that: (i) Xa < Xb (resp. Xa > x^); (ii) 
there is no critical line x = x^ such that Xa < Xc < x^ (resp. Xa > Xc > Xt). 
Also, let Pa be a real point of C lying on the line x = Xa, generated by ta G M, 
and let Vb = {^6,1, . . . , tfe,nb7 } (including —00, +00, if Poo belongs to the line 
X = Xb) be the set of real values generating the real points of C (1 {x = Xb}. 
The following statements are true: 

(1) If x'{ta) > 0, then Pa must be connected with the point Pb of C C] {x = Xb} 
generated by the least (resp. greatest) element ofVb which is greater (resp. 
less) than ta. 

(2) If x'ita) < , then Pa must be connected with the point Pb of C \~\ {x = Xb} 
generated by the greatest (resp. least) element of Vb which is less (resp. 
greater) than ta. 

Proof. We prove (1) for the case when Xa < Xb, the proofs of (1) for the case 
Xa > Xb, and of (2) in both cases, are analogous. Now let tc G Vb be the least 
element of Vb which is greater than ta- Since by hypothesis C has no vertical 
asymptotes, then Pa must be connected either with exactly one real point of 
Cr\{x = Xb} generated by a real value of the parameter, or with Poo- Now, we 
distinguish two different cases, depending on whether Poo belongs to the line 
X = Xb, or not. We begin with the case when Poo does not belong to x = Xb- So, 
Pa is connected with a point of C H {x = Xb} generated by some i E Vb- First 
of all, observe that t > ta- Indeed, by hypothesis C has no vertical asymptotes. 
Then, x{t) is defined for every t between ta and t, and since x{t) is a quotient 
of polynomials, x'{t) is also differentiable there. Moreover, x'{t) cannot vanish 
between ta and t because by hypothesis there does not exist any critical line 
between x = Xa and x = Xb- Hence, the sign of x'{t) is constant in the interval 
lying between ta and i, and since x'(ta) > 0, then x'{t) > in that interval; 
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therefore, x{t) is increasing there. So, since Xa = x{ta) < x{t) we deduce that 
ta < i. 

Now our aim is to prove that t = t^. For this purpose, observe that t > 
because t^. is the least element of Vh greater than ta'-, hence, we just have to 
prove that i > cannot occur. Assume by contradiction that i > tc- Since 
x{tc) = Xb = x{t) and x{t) is differentiable along [ta,t), by RoUe's Theorem 
x'{t) must vanish at some point of {tc, i). However, this is absurd because x'{t) 
is strictly positive in [ta,t), which contains {tc,t). 

Finally, let us consider the case when Poo belongs to the line x = Xb. If there 
exists t EVb, t +00, with t > ta, then P must be connected with P = (p{t), 
since otherwise by adapting the above argument one reaches a contradiction. 
On the other hand, if ta is greater than every real element of Vf,, then Pa 
cannot be connected with any other point of C fl {x = Xb} but Poo] however, 
since we consider Poo generated by +00, and ta < 00, the rule also holds in 
this case. □ 

Y 




Line 



Fig. 1. Connecting Points 



3.4 Full Algorithm 

The following algorithm Planar-Top can be derived from the preceding sub- 
sections. 

Planar-Top Algorithm: 
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Input: a planar curve C, parametrized by 



ip{t) = ixit),yit)) 



fulfilling: (i) Pi{t),qi{t) G Z[t] for i = 1,2, gcd{pi,qi) = 1 for z = 1,2; (ii) ip{t) 
is proper; (iii) C has no vertical asymptotes. 

Output: a planar graph Q homeomorphic to the curve. 

(1) (Critical Points) Compute the polynomial m{t) in Corollary [5l and the real 
roots of m(t). Then, compute: 

(1.1) The critical points of C (by evaluating x{t), y{t) at the real roots of m{t)). 
Store these points in a list 

For each of these points, store its coordinates, and the list of real t-values 
generating them. 

(1.2) The point Poo = (a^oo, Uoo) (if it exists), and the list of t-values generating 
it. 

(2) (Points of C on Critical Lines) 

(2.1) For I from 1 to r, compute the real points of C lying on the line x = x^. 
Store these points in a list 



£f — [Pt,!-, ■ ■ ■ ,Pl 



For each of these points, store its coordinates, and the list of real t-values 
generating them. 

(2.2) Check whether Poo belongs to some of the above lines x = xg. In the 
affirmative case, go to (3); otherwise, compute the real points of C lying 
on X = Xoo- Store these points in a list 

£oo [-foOjl; • • • ; -Poo,m] 

For each of these points, store its coordinates, and the list of real t-values 
generating them. 
(3) (Points of C on Non-Critical lines) 

(3.1) Let M = {ai...,as}, ai < ... < Og, be the set consisting of the x- 
coordinates of the critical points computed in (1.1) and (1.2). Also, let 
ao = Oi — 1, (is = tts + 1, and for i from 1 to s — 1 let aj = (aj + aj+i)/2. 

(3.2) For j from to s, compute the real points of C lying on the line x = cij] 
store these points in a list 



A/} = [P,,i,...,/^ 



For each of these points, store its coordinates, and the list of real t-values 
generating them. 
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(4) (Edges) 

(4.1) For i from to s — 1, connect the points of C lying on x = di and x = Oj+i 
by applying Theorem [SI 

(4.2) For i from 1 to s, connect the points of C lying on x = di and x = a, by 
applying Theorem [SI 

(5) (Isolated vertices) Compute the real isolated points of the curve, and add 
them to the graph. 

We will provide several examples of the output of the algorithm in Section \5\ 



4 Computation of the Graph Associated with a Space Curve 



In this section we let C C M'^ be a real curve, parametrized by 

'Plit) P2{t) P3{t)\ 



ifit) = {x{t),y{t),z{t)) 



.giWg2(t)'g3(t), 



where gcd(pi,gi) = gcd(p2,Q'2) = gcd(p3,g3) = 1. In the sequel, we consider 
the problem of algorithmically computing a graph Q homeomorphic to C. In 
order to do that, we will follow the strategy used to address the implicit case 
in (0)5 (0), (@); more precisely, we need to perform the following steps: 

(1) Project the curve onto a coordinate plane (the xy-plane, in our case) 

(2) Compute the graph Q associated with the projection (by using the algorithm 
given in Section [3]) 

(3) Lift the graph Q of the projection, to get Q. 

As in (0) and (0), here we will require just one projection in order to perform 
the lifting phase. Now in the following subsections we first describe the hy- 
potheses that we request on the input curve (essentially, that it is properly 
parametrized, and that it is "correctly placed" in space); then, we present the 
ideas and results needed for computing the vertices and edges of the graph, 
and finally we provide the full algorithm. 



4-1 Hypotheses 



Since C is rational, if it is not a line parallel to the z-axis, then its projection 
onto the xy plane, denoted as 7ixy{C), is an algebraic rational curve and can 
be parametrized by 



iP{t) = {x{t),y{t)) 



' Pi{t) pM 
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Thus, in the sequel we assume that the following hypotheses hold: 

(i) C has no asymptotes parallel to the 2;-axis (in particular, it is not normal 
to the xy-plane). 

(ii) ip(t) is a proper parametrization of n^y^C). 

(iii) TTxyiC) has no asymptotes parallel to the y-axis. 

In particular, hypotheses (ii) and (iii) imply that the graph of vr^y(C) can be 
computed by using the Planar-Top Algorithm. Now if the parametrization 
(p{t) of C is not proper, then ip{t) cannot be a proper parametrization of 
7Txy{C) either; then, in particular (ii) implies that C is properly parametrized. 
However, the converse does not necessarily hold, i.e. it can happen that (f{t) 
is proper, but ip{t) is not. From Section [3|, we know how to check hypotheses 
(ii) and (iii). In order to check hypothesis (i), the next lemma, analogous to 
Lemma O can be applied. 

Lemma 9 C has an asymptote parallel to the z-axis iff one of the following 
conditions happen: (a) qsit) has some real root, which is not a real root of 
Qiit) ■ q2{t); (h) deg(p3)) > deg(g3) hut deg(p2) < deg(g2), deg(pi) < deg(gi). 

Moreover, hypothesis (i) implies the following relationship between the points 
Qoo = liint^itoo V'(^)5 -Poo = linit^ioo V'(^)- Recall from Section [2] that they are 
the only points of HxyiC) and C, respectively, that may not be reached by any 
complex value of the parameter. 

Lemma 10 Assume that hypothesis (i) holds. Then, exists iff Qoo exists, 

and TXxy{Poo) = Qoo- 

Proof. If Poo exists, then it is clear that Qoo exists and is the projection of 
Poo- Conversely, if Qoo = (a^oo^^oo) £ then Poo exists because C has no 
asymptotes. □ 

On the other hand, hypothesis (ii) leads to the following result. Here, the 
notion of birationality arises; essentially, the projection of C is said to be bi- 
rational if there are not two different branches of C whose projections overlap 
(see Chapter 5 in ([^ for further information on birationality). 

Theorem 11 Assume that C is not a line parallel to the z-axis. Then, if ip{t) 
is proper, the projection of C onto the xy-plane is birational. Conversely, if 
(p{t) is proper and the projection of C onto the xy-plane is birational, then 
i/j{t) is proper. 

Proof. Let us see (^). For this purpose, let Q G iixy{C), Q ^ Q^o, satisfying 
that there are at least two different points P,P & C projecting onto Q. Since 
Q 7^ Qoo, by Lemma [10] none of these points is Poo, and hence both are 
reached by (p{t). Let tp ^ tp be the t- values generating P, P, respectively. Then, 



13 



ip{tp) = ip(tp), and thus Q is generated by two different values of the parameter. 
But since ip{t) is proper, this can only happen for finitely many points, and 
thus the projection is birational. Conversely, given any Q e n^y^C), Q ^ Qoo, 
not generated by any root of qsit) (notice that we are excluding finitely many 
points), the t- values reaching Q are exactly those ones generating the points 
of C that are projected onto Q. Since (p(t) is proper, almost all points of C 
are generated by just one value of the parameter. And since the projection is 
birational, we conclude that almost all points of iixyiC) come from just one 
point of C, and therefore almost all points of TixyiC) are generated by just one 
value of the parameter. So, (<^=) holds. □ 

It is well-known that almost all affine transformations of the type {X = x + 
az, y = y + bz, z} transform C so that its xy-projection is birational. So, if ip{t) 
is proper, almost all of these transformations set C proper. Moreover, if ip{t) is 
not proper there exist reparametrization algorithms (see (jsj), (IigI)). Therefore, 
in the sequel we will assume that the above hypotheses hold. 

4-2 Definition of the Space Graph. 

By assuming the hypotheses of the preceding subsection hold, we can compute 
the graph Q associated with iTxyiC) with the Planar- Top Algorithm described 
in Section [31 Hence, in the following we will assume that this process has 
already been carried out. 

Now we make precise the definition of the graph Q that we want to compute. 

Definition 12 LetC be a space curve in the above conditions. Then, the graph 
associated with C, Q , is the following graph: 

(i) Its vertices are the real points of C giving rise (by projection) to the vertices 
ofG. 

(a) Its edges are the result of "lifting" to space the edges of Q , i.e. of computing, 
for each edge i of Q, an space segment i' corresponding to the branch of C 
giving rise (by projection) to i. 

Hence, we have to lift to space the vertices and edges of Q in order to compute 
Q. Let us see that this lifting operation is well-defined. 

Theorem 13 Every vertex of Q , except perhaps the isolated vertices, lifts to 
at least one real space point of C. 

Proof. Every point of Q G T^xyiC) HM^ fulfills one of the following conditions: 
(1) Q = Qoo', (2) there exists to satisfying that Q = ip(to); (3) Q does not 
fulfill (2), but there exists to & C such that Q = ipito)- the first ca.se,Q is 
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lifted to Poo by Lemma [TOl If Q belongs to the second group, then it is lifted 
to P = <y9(to) because C has no asymptotes parallel to the 2;-axis. Finally, if Q 
belongs to the third group then it is an isolated point of HxyiC)] in this case, 
Q comes from a real point of C iff 2;(to) G IR- □ 

Remark 2 Real isolated points of 7ixy{C) may come from real isolated points 
of C, or from points of C whose z-coordinate is complex. In any case, thanks 
to hypothesis (i) they do not come from branches of C normal to the xy -plane. 

Now let us consider the lifting of the edges of the planar graph. The next result 
guarantees that, under the considered hypotheses, this lifting process can be 
always carried out. In particular, it implies that there are no real branches of 
T^xyiC) coming from complex components of C (which is a phenomenon that in 
general can happen when working with space algebraic curves; see for example 
p. 734 in (0)). 

Theorem 14 Under the considered hypothesis, for every edge I of Q there 
exists one and just one branch of C giving rise to i. 

Proof. Let £ be an edge of Q. By construction of the graph provided in Section 
m , if Qoo exists, it is always included as a vertex of ^. So there exists a real 
open interval J C M such that generates the real branch of TTxy{C), that 
we denote by £, corresponding to I. On the other hand, for every t ^ I we 
have that z{t) must be defined, because otherwise C has an asymptote parallel 
to the z-axis. Then, ip{t) is defined for every t & I, and gives rise to a real 
connected branch of C projecting as C Furthermore, since ip{t) is proper by 
hypothesis, the projection onto the xy-p\a.ne is birational by Theorem [TTl 
Hence, there are just finitely many points of C giving rise, by projection, to 
the same point of iTxyiC); but none of these points can give rise to a point of 
£, because such a point would create a singularity of HxyiC) which would split 
£ into two different edges, and £ is already an edge of Q. Then, we conclude 
that C lifts to a unique connected real branch of C. □ 

4.3 Computation of the Vertices 

From Definition [121 this process is the lifting of the vertices of Q. From the 
construction of the planar graph, one may see that for each vertex Qi = 
{xi,yi) of Q the algorithm stores the real values tj^i, . . . ,ti^r of the parameter 
generating it. For a fixed i, z{tij) is well-defined for j e {l,...,r}, since 
otherwise C has an asymptote parallel to the z-axis. Hence, Qi is lifted to the 
space points 

Furthermore, if Qoo exists, then it is hfted to Poo and to the space points 
reached by the real values of the parameter generating Qoo, if any. Proceeding 
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this way, the only remaining space vertices are the real isolated ones (which, 
by Proposition [21 are generated by complex values of the parameter). So, in 
the rest of the subsection we consider this kind of points. 

From Theorem [131 the real isolated vertices of vr^y(C) do not necessarily come 
from real isolated points of C (since they may be the projection of complex 
space points). Conversely, a real isolated point of C does not necessarily project 
as an isolated point of nxy{C), because its projection may coincide with the 
projection of some other real point of C which is not isolated. However, the 
next result ensures that isolated points of C always project as vertices of Q; 
therefore, these points are computed when lifting the planar vertices. 

Lemma 15 Let P E C he a real isolated point. Then, 7ixy{P) is a vertex of Q . 

Proof. If Tx,j.y{P) is an isolated point of 7r^.y(C), then the statement is true. 
Otherwise, there exists a point P' 7^ P in a real branch of C such that 
TTxyiP) — '^xy{P')- Observe that P cannot be Pqo because it is isolated. There- 
fore, suppose that it is reached via ip{t) by tp G C. Now we distinguish the 
cases P' 7^ Poo or P' = Poo, respectively. If P' 7^ Poo, then P' = <y5(V) ^'^^^ 
tp/ G M. Thus, Tixy{P) is generated via ip{t) by two different values of the pa- 
rameter, namely tp, tp/, and since ip{t) is proper, iixyiP) is a self-intersection of 
TixyiC). Hence, it is a singularity of 'KxyiC), and the statement follows. Finally, 
if P' = Poo then 7ixy{P') = Qoo and therefore it is also a vertex of Q. □ 

Then, we might recover isolated singularities of C by determining the complex 
values of the parameter that generate (by projection) vertices of and by 
computing those real points of C which are generated by those values. Nev- 
ertheless, in the sequel we consider an alternative method, analogous to that 
in Subsection 13. 2[ For this purpose, the following lemma is needed. Here, we 
denote a complex value of the parameter t as t = u + iv., where = — 1 and 
M, w G M. Also, we write 



and 



Piiu + iv) ■ qi{u-\- iv) 


1 






P2{u + iv) ■ q2{u + iv) 


1 






Ps{u + iv) ■ q3{u + iv) 


1 







(a(M, v) + ib{u, v)) 
(c(m, v) + id{u, v)) , 

{e{u, v) + ih{u, v)) 



Then the following result, analogous to Lemma [3, holds. Here, q{u, v) denotes 
the result of substituting t = u + iv in lcm(gi, Q's)- As in Lemma U\ by 
applying the following result one computes a finite set of complex points which 
contains the complex points generating the isolated singularities of the space 
curve. 
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Lemma 16 Let P G C fl M^. Then, P is generated by a complex value of 
the parameter t^ = uq + Wq if and only if there exists G K satisfying that 
{uo,Vo,Wo) is a real solution of the system 



b{u,v) = 
d(u, v) = 

(3) 

h{u, v) = 
V ■ \q{u, f )p ■ ii; — 1 = 

As in the planar case, one can certify the number of real solutions of the system 
by Hermite's method; also, one can construct another system whose solutions 
correspond to complex values of the parameter generating points that are also 
reached by real values of the parameter, and proceed as in the 2D case. 



4-4 Computation of the Edges 



The method consists of the lifting of the edges of Q. So, let £ be an edge of Q; 
by Theorem UM. i is lifted to an space edge £' G Q. In order to compute the 
crucial observation is that the computation of the edges of Q is in fact done 
by connecting not points, but values of the parameter t. Hence, each edge £ 
can be identified with a pair 

where ta,t belong to M U { — oo, +oo}, and where the vertices of Q defining £ 
are Qa = '^{ta), Q = i^it) (see also Figure 1); here, Qoo = '?/'(±oo). Hence, £ is 
lifted to the space segment connecting the points Pa = fita), P = also, 
-Poo = (^(ioo). Notice that this idea works perfectly when Q is the projection 
of several real points of C. 



4-5 Full Algorithm 



The following algorithm Space-Top can be derived from the preceding sub- 
sections. 



Space-Top Algorithm: 



Input: a space curve C, parametrized by 

'pi{t) P2{t) pM 



^{t) = {x{t),y{t),z{t)) 



Q2{ty q2{ty qz{t), 
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fulfilling: (i) q,{t) G Z[t] for 2 = 1,2, 3, gcd(p„ g^) = 1 for z = 1, 2, 3; (ii) C 

' pijt) P2{ty 



has no asymptotes parallel to the z-axis; (iii) ■i}j{t) = ( ^-tt- 1 is proper; 



(iv) T^xy{C) has no asymptotes parallel to the y-axis. 
Output: a space graph Q homeomorphic to C. 

(1) (Projection) Compute the graph Q associated with the projection nxy{C) of 
C onto the xy-plane, parametrized by ipit), by applying Planar-Top. 

(2) (Lifting phase) 

(2.1) (Vertices) For each vertex of Q: if P is generated by ti,...,tp where 
ti, . . . ,tp G M U {—00, +00}, then P lifts to the points (p(ti), . . . , f(tp); 
Qoo, if it exists, lifts to Poo, and we write Pqo = v^(±C)o). 

(2.2) (Edges) For each edge of Q: if £ is identified (according to Subsection 14.41) 
with [ta,tb], where ta,tb G M U {—00, +00}, then it is hfted to the space 
edge obtained by connecting (p(ta), ^{h) by means of a segment. 

(3) (Isolated vertices) Add to Q the real isolated singularities of C. 

Several examples of the output of this algorithm are presented in the next 
section. 



5 Experimentation and Examples 



The algorithm has been implemented in Maple 13, and the examples run 
on an Intel Core 2 Duo processor with speeds revving up to 1.83 GHz. The 
implementation allows the option of computing isolated points, or not. The 
reason for introducing this option is that the number of isolated points is 
certified by means of Hermite's method, and this method may be costly. 

On the other hand, the user can decide the number of digits used in the 
computation. Suppose we denote such a number by n. Then, when running 
the algorithm, the computing starts using n digits. However, if the algorithm 
detects that the number of points in a vertical line is not the right number, the 
precision is automatically increased by 5 more digits and the whole process 
starts again. In our experimentations, we usually set n = 10, the default value 
of Digits variable in Maple, and in the implementation, the number of digits 
is limited to a maximum of 500, although we have never needed more than 70 
digits. 

Next, we first present examples of the 2D algorithm. In Table 1, we include, 
for each curve, the degree of the parametrization (i.e. the maximum exponent 
of the parameter in the numerators and denominators of the components of 
the parametrization, dp), the total degree of the implicit equation (dt), the 



18 



number of terms of the implicit equation (n. terms), the timings in seconds 
corresponding to the graph without computing isolated points (to) or com- 
puting them (ti), and the number of digits used in the computations. The 
parametrizations corresponding to these examples are given in Appendix I. 



Example 






n. terms 


tn 


^ i 


Digits 


1 


Q 

o 


u 


1 u 


yj.oocf 




1 n 


9 


Q 

o 


fi 
o 






1 078 


1 n 




8 


8 


q 




71 1 72 


40 


4 


4 


4 


7 


0.109 


0.110 


10 


5 


6 


6 


28 


0.203 


11.859 


10 


6 


8 


8 


21 


0.171 


2.50 


10 


7 


23 


23 


335 


49.797 


> 1 h. 


10 


8 


6 


12 


49 


13.625 


> 1 h. 


10 


9 


17 


17 


171 


1.656 


> 1 h. 


10 



Table 1: 2D Examples. 



One may notice that as the degree increases (see Example 7 or Example 9) the 
computation of the isolated points turns very costly. An alternative for those 
cases could be to detect isolated points directly by checking the existence of 
complex values of the parameter corresponding to real singular points; users 
interested in certifying rigourously the number of isolated singularities, can 
choose to apply Hermite's method later. 

The pictures corresponding to the examples in Table 1 can be found in Figure 
2; from left to right we have Examples 1, 2, 3 in the first row, 4, 5, 6 in the 
second row and 7, 8, 9 in the third one. Examples 2 and 6 are the offsets 
of the cardioid and of the cubical cusp, respectively; furthermore. Example 4 
corresponds to the epitrochoid. Notice that the curves in Examples 2, 3 and 
4 are not in generic position. 

Finally, we present examples of the 3D algorithm. In Table 2, for each curve we 
include: the degree of the parametrization {dp), the total degree of the implicit 
equation of the projection onto the x|/-plane (cij), the number of terms of this 
projection (n. terms), the timing without computing isolated points (to), the 
timing including the computation of isolated points (ti), and the number of 
digits used. As in the 2D-case, in all cases the computations start with 10 
digits, and the algorithm increases the number of digits when it is necessary. 
The parametrizations corresponding to these examples are given in Appendix 
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Fig. 2. Examples of the 2D algorithm. 



II. 
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Example 


dp 


di 


n. terms 


to 


ti 


Digits 


1 


8 


8 


38 


5.578 


6.188 


30 


2 


10 


10 


65 


3.516 


3.297 


10 


3 


21 


21 


234 


4.453 


4.515 


10 


4 


4 


7 


8 


0.657 


0.625 


10 


5 


6 


6 


28 


0.437 


0.266 


10 


6 


8 


4 


5 


0.141 


0.109 


10 


7 


4 


4 


15 


0.125 


0.500 


10 


8 


12 


12 


91 


1.015 


0.875 


10 


9 


16 


16 


142 


74.00 


74.094 


65 



Table 2: 3D Examples. 



The pictures corresponding to these curves can be found in Figure 3. The 
diamond in each picture points out the origin of the system of coordinates; 
moreover, in Example 7 we have not plotted the axes for the isolated point to 
be better appreciated. 
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Fig. 3. Examples of the 3D algorithm. 



References 



[1] Arnon D., MacCallum S. (1988). A polynomial time algorithm for the 
topology type of a real algebraic curve, Journal of Symbolic Computation, 
vol. 5 pp 213-236. 

[2] Alcazar J.G., Sendra R. (2005) Computation of the Topology of Real Al- 
gebraic Space Curves, Journal of Symbolic Computation 39, pp. 719-744. 

[3] Andradas C, Recio T., Sendra J.R. (1997) A relatively optimal 
reparametrization algorithm through canonical divisors. Proceedings 
ISAAC 97, ACM press, pp. 349-355. 

[4] Andradas C, Recio T. (2007) Missing points and branches of real para- 



22 



metric curves, Applicable Algebra in Engineering and Computing 18 (1- 
2), pp. 107-126 

[5] Cox D., Little J., O'Sliea D. (1992). Ideals, Varieties and Algorithms. 
Springer. 

[6] Cox D., Little J., O'Sliea D. (2005). Using Algebraic Geometry. Second 
Edition. Springer. 

[7] Diatta D.N., Mourrain B. and Ruatta O. (2008) On the Computation 

of the Topology of a Non-Reduced Implicit Space Curve. In Proceedings 

ISSAC 2008, ed. David Jeffrey, pp. 47-55. 
[8] Eigenwilling A., Kerber M., Wolpert N. (2007) Fast and Exact Ceometric 

Analysis of Real Algebraic Plane Curves, in C.W. Brown, editor, Proc. 

Int. Symp. Symbolic and Algebraic Computation, pp. 151-158, Waterloo, 

Canada. ACM. 

[9] El Kahoui M. (2008) Topology of Real Algebraic Space Curves. Journal of 

Symbolic Computation vol. 43, pp. 235-258. 
[10] Gianni P., Traverso C. (1983). Shape determination of real curves and 

surfaces, Ann. Univ. Ferrera Sez VII Sec. Math. XXIX pp 87-109. 
[11] Gonzalez- Vega L., El Kahoui M. (1996). An improved upper complexity 

bound for the topology computation of a real algebraic plane curve, J. 

Complexity 12 pp 527-544. 
[12] Gonzalez- Vega L., Necula I. (2002). Efficient topology determination of 

implicitly defined algebraic plane curves. Computer Aided Geometric De- 
sign, vol. 19 pp. 719-743. 
[13] Hong H. (1996). An effective method for analyzing the topology of plane 

real algebraic curves. Math. Comput. Simulation 42 pp. 571-582 
[14] Perez-Diaz S. (2007). Computation of the singularities of parametric plane 

curves. Journal of Symbolic Computation, vol. 42, pp. 835-857. 
[15] Rubio R., Serradilla J.M., Velez M.P. (2008). Detecting real singularities 

of a space curve from a real rational parametrization. Journal of Symbolic 

Computation, etc. 

[16] Sederberg T.W. (1986). Improperly parametrized rational curves. Com- 
puter Aided Geometric Design 3, 67-75. 

[17] Seidel R., Wolpert N. (2005) On the Exact Computation of the Topology 
of Real Algebraic Curves. Proc. of the 21st Ann. ACM Symp. on Comp. 
Geom. (SCG 2005). ACM, 2005 107-115. 

[18] Sendra J. R. (2002). Normal Parametrizations of Algebraic Plane Curves. 
Journal of Symbolic Computation vol. 33, pp. 863-885. 

[19] Sendra J.R., Winkler P., Perez-Diaz P. (2008). Rational Algebraic Curves, 
Springer- Verlag. 



23 



Appendix I: Parametrizations of the planar curves used in the ex- 
perimentation 



Example 1: 

'44 + 37t^ - 23^2 + 87t 95 - 61t^ - St^ - 29t 



^ 10 + 29t^ + 98t2 - 23t ' 40 + - 49^^ - A7t J 
Example 2: 

_ / 3456^5 - 31104^3 + 6f - 756f + 61236^2 - 39366 
~ V 486^4 + 36^6 + 2916^2 + t^ + 6561 ' 

-18t(864t3 - I6t^ - 1296t + 6^^ - 126t^ - 1134t2 + 4374) \ 
486t4 + 36t6 + 2916^2 + t^ + 6561 ) 

Example 3: 

^{t) = {t^ - 8t^ + 20t^ - 16^2 + 2,f - 7t^ + Ut^ - 7t) 
Example 4: 
, , / -7t^ + 288t2 + 256 -m^ + 256t \ 



^ + 32t2 + 256 ' + 32t2 + 256^ 
Example 5: 

' 3t2 + 3t + l t^{t^-2t + 2) 



^-3t - 2t4' -st-l+t^ - 2t^ J 

Example 6: 

, , / (^2 - l)(t4 _ 1 + 9^2) _(^8 _ + 2t2 - 1 - 54^4) 

= - 



9t2(t2 + l) ' 27(^2 + 1)^3 J 

Example 7: 

cp{t) = (-83^23 + 98^20 _ 48^18 - 19ti3 + Q2e^ + 37^8, _13 _ 54^2^ + 64t25 _ 90t22 _ QQfi^ 
Example 8: 

' 9 + 85t^ + 80t^ + + 74^2 + 27t -56 - 5t^ + 36t^ - 8t^ + 30^^ - 3t \ 



^5 - 91t6 + 81t^ + 65^4 - 12^2 + 78t' -79 - 70^^ + 42^^ + 9^3 - 21^2 - 27t J 
Example 9: 

(^(t) = + 80 - 20t^ - At^ - 89t^ - 77^2 + 69t, - 64 - 33^^ + 21^"^ - 35^^ + 97^2 + 30t) 
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Appendix II: Parametrizations of the space curves used in the ex- 
perimentation 



For each example, will use the notation (p{t) = {x{t),y{t), z{t)). 
Example 1: 

, , 36t(-l - 98t - 3954^2 - 78868^^ - 726692t^ - 1092840t^ + 31242296^*^ + 193263952t^ 
xit) = 



yit) 

z{t) 



lit) 

-648^2(1 + 84t + 2940^2 + 54880t^ + 549996t^ + 2492112^^ + 2385712^^) 

36t(476t3 + 426^2 + 42t + 1) 
64660^4 + 10976^3 + 1176^2 + set + 1 ' 



with = 112t + 5488^2 + 153664^^ + 2741608t^ + 33057472t5 + 272552896*^ + 
1419416320t^ + 4180915600t*^ + 1 

Example 2: 

, , 7 + 33^1° + 80^9 - 57 f + 88t^ + 75^^ 
x{t) 

y{t) 

z{t) 



5^10 + 61t8 + 8f + 71t6 - 16^5 + 37t 

18^8 + 28t^ + 58t^ + 69t^ + 8^^ + 4t 
5tio + 61^8 + 8f + 71t^ - 16t5 + 37t 

-94^9 - 59t^ + - 82t^ + 69^2 - t 
5tio + 61t8 + 8f + 71t^ - 16t^ + 37t 

Example 3: 

'^^^^ "1^ + 1 ' t2 + 1 ' t2 + 1 j 

Example 4: 
Example 5: 

(t- 1)4(1 + 4t + 7t2) 



xit) 

yit) 



1 - 4t + 17^2 - 5t6 - 13t4 + 20^5 + 48^3 

(1 - 4t + 22^2 - At^ + t%l + tf 
1 - 4t + 17t2 - 5t6 - 13^4 + 20^5 + 48^3 
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_ (l-4t + 22t^-4t3 + t^)(l+t)^ 
^ ~ 1 - 4t + 17^2 - 5t6 - 13^4 + 20t5 + 48t3 

Example 6: 

- 1 - (1 - 3t^)t (l-3t2)tn 

Example 7: 

57 - 7t^ + 22^3 - 55^2 - 94t 



y{t) 

z{t) 



-73 - 56^4 - 62t2 + 97t 

-82 - - 83t^ - 10^2 + Q2t 
-73 - 56t^ - 62^2 + 97t 

-82 - 4^4 - 83t3 - 10^2 + Q2t 



-73 - 56^4 - 62^2 + 97t 
Example 8: 

x{t) = 91 + 11^12 _ 49^10 _ 47^7 ^ 4Q^6 _ 
y{t) = -28t^2 ^_ iQ^io ^_ 30^8 _ 27^5 - 15t^ - 59t'^ 
z{t) = 53 + 43^1° + 92^9 - 91^6 - 88^^ - 48t 
Example 9: 

x{t) = -m^^ + 81t^ + 65t^ - I2t^ + 78t^ + 5t^ 
y(t) = -70t^^ + 42ti5 + 9ti2 _ 21^9 - 27t« - 79f^ 
z{t) = 62 - 14^14 + 83^12 - 96^7 _ st^ _ 54^2 
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